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A unitary framework based on the Bakamjian- Thomas 
construction of relativistic quantum mechanics is used to de- 
scribe two-pion scattering from threshold to 1400 MeV. The 
framework properly includes unitarity cuts for one-, two- and 
three-hadron states and provides an excellent description of 
the available data for nn phase shifts and inelasticities. The 
role and importance of three-hadron cuts are calculated and 
discussed. 



I. INTRODUCTION 

A nonperturbative framework capable of describing the 
relativistic, coupled-channel scattering of hadrons is pre- 
sented. The approach is based on a relativistic Hamil- 
tonian formulation with model interactions introduced 
into the mass operator, and with few-body states imple- 
mented in a way that maintains the unitarity of the the- 
ory. The elementary degrees of freedom in the framework 
are finite-sized hadrons which provide a natural ultravi- 
olet regularization, ensuring that the scattering ampli- 
tudes are finite. 

The Hilbert space is truncated to include only one-, 
two- and three-body states. A central and novel fea- 
ture of this framework is the explicit inclusion of both 
real and imaginary parts of scattering amplitudes arising 
from the opening of three-body channels. The proper 
handling of three-body unitarity cuts is crucial to gain- 
ing a deeper understanding of several well-known scatter- 
ing systems; a good example is nN scattering in the P\\ 
channel, which exhibits a significant inelasticity arising 
from the intermediate three-body irirN state . 

It presents a formidable challenge to develop a general, 
relativistic scattering framework to describe the final- 
state interactions between hadrons, that includes the ef- 
fects of three-body unitarity cuts. Nonetheless, a practi- 
cal framework which can treat hadron reactions beyond 
the lowest-order valence quark picture is clearly desirable. 
For example, the systematic analysis of hadron reactions 
in the baryon resonance region currently being conducted 
in Hall B at the Thomas Jefferson National Accelerator 
Facility (TJNAF) requires that such a framework be used 
to extract information about baryon resonances in this 
highly complex dynamical region. The framework devel- 
oped herein is an attempt to construct a useful, relativis- 
tic framework capable of describing the nonperturbative, 
low-momentum transfer final-state interactions between 
hadrons in a unitary manner. 

For the first application of the framework developed 
here, a simple model for 7T7t scattering is introduced and 
used to described the S and P partial waves for energies 



ranging from the two-pion threshold up to 1400 MeV. 
This system provides an excellent test for the framework. 
A relativistic treatment is quite important when dealing 
with particles as light as pions and the interplay between 
strong dynamics and chiral symmetry makes this system 
quite interesting. The 7T7t system is somewhat simpler 
than others, in that its study requires only a minimal 
complication from the proper implementation of rela- 
tivistic spins, since both of the two-body states (7T7T and 
KK) involved are comprised of spin-0 particles. Another 
attractive aspect of applying the framework to 7T7t scat- 
tering is the relative wealth of experimental data for the 
isoscalar, S- wave channel. 

One drawback with using 7T7t scattering as a touchstone 
may be that inelasticities due to open states of three or 
more particles, do not appear to be significant for this 
process; that is, the pnir and 77777777 thresholds seem to 
have little impact on the S- and P-wave observables. The 
effects of the opening of these three- and four-body chan- 
nels seem to be overwhelmed by the opening of the two- 
body KK channel. Nonetheless, several important as- 
pects of the framework can be explored in an application 
to 7777 scattering. 

The isoscalar-scalar (1 = 0, J PC = ++ ) channel of 7T7T 
scattering has been a subject of numerous and extensive 
studies. The study of meson scattering in this low-energy 
region may be an ideal test of our understanding of the 
interplay between bound states in QCD and chiral dy- 
namics. The region near E = 1000 MeV is perhaps most 
interesting, as it is dominated by the mixing between 7777 
and KK channels and the isoscalar-scalar /o (980) meson 
resonance. The nature of the /o(980) resonance, and the 
question of whether it is comprised of valence quarks or 
arises purely from meson scattering dynamics, has been 
addressed by many authors [^|-^| . Above this energy re- 
gion, three additional scalar meson resonances have been 
well established. These are referred to as the /o(1370), 
/o(1500) and the /o(1710). It is still unclear which of 
these should be considered as quark-antiquark bound 
states, glueballs or possibly resonances arising from dy- 
namical effects of final-state interactions p|-p^]. Thus, 
the isoscalar-scalar channel remains a source of great in- 
terest and mystery for meson phenomenology. 

Although there have been previous studies which em- 
ploy a framework similar to the one developed here, there 
are some important differences. Most studies of meson 
scattering dynamics are based on potential models (as is 
the framework developed here.) However, most other ap- 
proaches typically include one- and two-particle channels 
only; that is, they include s-channel states and several 
two-particle channels, such as 7777, KK, o~a, etc. They 
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either neglect the possibility of open three-particle chan- 
nels altogether or only partially implement them. For ex- 
ample, in the model developed by the Julich group [ O , 
the interaction potentials between two-particle channels, 
such as tttt-tttt or tttt-KK interactions are obtained us- 
ing an instantaneous approximation of a meson-exchange 
model. Such instantaneous approximations generally do 
not account for absorptive effects due to the opening of 
three-body channels. Still, there is no question that the 
Julich model quite successfully describes the phase shifts 
and inelasticities of tttt scattering for S, P and D waves. 
Alternatively, the Krakow group jl5| has developed a 
separable-potential model for tttt scattering. In their cal- 
culation, few-body dynamical effects are incorporated by 
including additional, effective two-body channels, such 
as a a a channel . Their model also obtains excellent 
results for the tttt phase shifts and inelasticities. 

The outline of this article is as follows. In Sec. 0, 
the relativistic scattering formalism employed herein is 
briefly discussed, beginning with a short proof of the co- 
variance of observables calculated within this framework. 
Then, the integral equations that relate the one-, two- 
and three-body scattering T-matrices are provided. 
In Sec. Ill, the framework is applied to a study of tttt 
The particle states that are included in the 
along with the necessary dynam- 
The interactions employed in 



scattering. 

model are discussed, 
ical model parameters, 
this study arise from the meson exchanges which cou- 
ple states of various numbers of particles to each other. 
In our framework, these interactions arise from one-, two- 
and three-meson intermediate states which may exhibit 
production thresholds, resulting in absorptive contribu- 
tions to the kernels and self-energies appearing in the 
Lippmann-Schwinger equations. A simple fitting proce- 
dure is shown to provide excellent agreement with data 
for tttt scattering phase shifts and inelasticities. Details 
of the relevant model dynamics that produce the various 
features observed in the resulting phase shifts and inelas- 
ticities are discussed. Then, it is shown that the numer- 
ical methods employed herein are sufficient to maintain 
the unitarity of the framework to better than one part in 
a million. Finally in Sec. IV, the article is summarized 



and plans for future studies are presented. 



II. RELATIVISTIC QUANTUM MECHANICAL 
FRAMEWORK 

In this section, a relativistic Hamiltonian framework 
that provides a covariant unitary approach to the study 
of multichannel scattering is described. Lorentz sym- 
metry is maintained by identifying the interactions with 
the mass operator (that is, the Hamiltonian in the over- 
all center-of-momentum frame). It is shown in Ref. p"^ ] 
that the complete set of Poincare generators can be con- 
structed in a simple way that separates the internal dy- 
namics from the center-of-momentum (CM) motion. In 



Sec. [I A a proof of the covariance of this approach is 
provided. Furthermore, Betz and Coester show that 
such a framework can satisfy cluster separability. All of 
these features are desirable for the study of hadron scat- 
tering. 
In Sec. 



HA, the Lorentz covariance of the framework 



is demonstrated and the fully-interacting mass operator 
M. is constructed. It is shown that the framework leads 
to Lorentz-invariant on-shell T-matrix elements T(_E, P) 
for colliding particles with total momentum P and en- 
ergy E] that is, one finds T(E,P) = T(s/s) where 
-y/i = \lE 2 — P 2 is the invariant mass of the system. 

The Poincare generators act on a Hilbert space which, 
in general, contains an infinite number of states. The 
Hilbert space is truncated to include only those states es- 
sential to describe the scattering system of interest within 
a particular energy range. Here, only one-, two- and 
three-particle states are maintained. Following this trun- 
cation, the operator form of the Lippmann-Schwinger 
equation can be written as a set of coupled integral equa- 
tions. The input that determines the dynamics is given 
in terms of the matrix elements of a model potential V. 
Once these are provided, the full scattering problem is 
solved in a straightforward manner. 



A. Relativistic covariance 

A simple realization of the Poincare algebra for an 
interacting system of a finite number of constituents 
is given by the Bakamjian-Thomas construction [fj"7| . 
This approach has the advantage of providing a Lorentz- 
covariant generalization for a large class of noncovari- 
ant microscopic models, such as the constituent quark 
model. In principle, a noncovariant microscopic model 
could be used to obtain matrix elements of the under- 
lying elementary hadronic potentials V. In this case, 
one might consider this framework as a means to extend 
the original noncovariant model dynamics, allowing for a 
Lorentz-covariant treatment of scattering phenomena. 

The explicit construction of the Poincare algebra pro- 
ceeds as follows. Starting from a system of noninteracting 
particles, described by their coordinates x a , momenta p a , 
spins s a , and masses m a , the Poincare generators are 



H = ^£(m Q ,p a ) ^^^ml + pl, 

a a 
a 

J = x p a + s a , 

a 



So X p a 



£ (m a ,p a ) + m a 



(2.1) 



Here, H and P are the total free energy and linear mo- 
mentum of the system, J and K are the total angu- 
lar momentum and boost operators, respectively. The 
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relative coordinates r Q , relative momenta k a , center-of- 
momentum (CM) spins s' a , and the CM coordinates R cm , 
total momentum P cm , and total spin S cm , are introduced 
via the Gartenhaus-Schwartz transformation which al- 
lows a separation of the internal dynamics and CM mo- 
tion. In terms of these new variables, the Poincare gen- 
erators are given by 



= v/P 2 +A* 2 (ki,k 2) ...), 
P — P 

* 1 cm: 

J — Pcm X Pcm S cm , 

1 S x P 

with the constraints, 

^ m a r a = 0, 

a 



(2.2) 



x k„ 



(2.3) 



In Eq. fl2.2|), the quantity M = M(ki, k 2 , . . .) is referred 
to as the free invariant mass in the Schrodinger picture. 
The internal momenta k a are related to the individual 
particle momenta p a via a free Lorentz transformation 
to the CM frame, 



k a = A(k a <- p a )Pa 
Pa P 



Pa 



M(M + H) 



£{m ai p a ) 
M 



(2.4) 



and the CM frame spins are related to the individ- 
ual spins s a via a Wigner rotation corresponding to the 
product of Lorentz boosts R = A(0 <— p a )A(p a <— 
k a )A(k a <- 0), leading to 



s' a = D^(R)s a D^(Ry 



(2.5) 



Interactions are incorporated into the Poincare gener- 
ators by the addition of a term in the free mass operator, 



M ^ A4/(r a ,k a ,s' a ) = M + V. 



(2.6) 



Thus, transforming the free Hamiltonian H into the in- 
teracting Hamiltonian Hi, 



H -> Hi = H + W, 

w = JmI + p 2 - Vm 2 + p 2 . 



(2.7) 



This replacement preserves the canonical commutation 
relations, provided V = V^(r a ,k a ,s' a ) is a function of 
internal coordinates only and is invariant under rotations 
[V,J] = [V, S cm ] = 0. For example, consider the case 
for which the elementary interaction is a Yukawa-type 



three-meson vertex. The matrix elements of the three- 
meson interaction vertex would be given by (a|V|6c) and 
would only depend on the internal variables associated 
with the CM frame where p a = p + p c = 0. Of course, 
these internal variables can be expressed in terms of the 
individual particle momenta in anoth er fr ame by u sing 
the boost relations analogous to Eqs. (2.4) and (2.5). 

Within this framework, the Lorentz covariance of ob- 
servables may be demonstrated from the following con- 
siderations. Construct an invariant T-matrix which sat- 
isfies a Lorentz-invariant Lippmann-Schwinger equation 
(LSE), 



r = v + vgr, 

with an invariant interaction V, 



(2. 



V = H 2 t — H 2 



(P 2 

W 2 
V 2 



<rM\) - (P 2 +M 2 ) 
- HW + WH 
MV + VM, 



(2.9) 



and an invariant propagator Q for scattering energy E, 
given by 



g = (E 2 - H 2 + it)- 1 

= {E 2 - P 2 - M 2 + ie)- 1 . 



(2.10) 



Since V is independent of the CM momentum P and the 
scattering energy E, one may rewrite the scattering en- 
ergy E = y/s + P 2 , in terms of a new variable s, referred 
to as the invariant mass squared. From Eq. ( |2.10 ), one 
observes that the propagator Q(s) is a function of the 
invari ant mass squared only, and one concludes that the 
LSE (2.8) depends only on the invariant mass squared 
s. It follows that the resulting T-matrix, T(y / s) depends 
only on the invariant mass squared s. 

It is possible to relate the on-shell matrix elements of 
this invariant T-matrix to the on-shell matrix elements 
of a T-matrix that is the solution of a non-invariant LSE 
with the interaction potential W, 

T{E,P) = W + WG(E) T(E,P), (2.11) 

where G(E) = (E — H -Me)" 1 . The relation between the 
on-shell matrix elements is given by 



T(x/i) = 2v/s + P 2 T(E,P), 



(2.12) 



which can be demonstrated term by term by expanding 
the on-shell matrix elements of T in powers of the po- 
tential V, 



T = V + VQ(E)V + 0(V 3 ) 



= W 2 + 2EW 
+W(E + H) 



1 



= 2EW + 2EW 
= 2ET(E,P). 



E 2 ~H 2 + n 

1 



-{E + H)W + Q{V 3 ) 



E-H+ie 



W + Q{V 6 ) 
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In this article, calculations are carried out in the CM 
frame for which E = y/s, then the interaction potential 
W = Mi — M = V, and the relevant LSE is 



T{E) = V + VG(E)T{E), 



(2.13) 



where T(E) = T(E, P = 0). In the CM frame, one finds 

T{y/a) = 2^~sT{yfs). (2.14) 

Thus, the on-shell matrix elements of the s olution T(E = 
y/s) of th e non-invariant LSE in Eq. ( 2.13| ) are related by 
Eq. (2.14) to the on-shell matrix elements of the solution 
T(yV) of the invariant LSE of Eg. ( [2§| ). It follows that 
observables calculated from Eq. (|2.13| ) are equivalent to 
those calculated from a Lorentz-invariant theory. 



B. Coupled Lippmann-Schwinger Equations 

In the above framework, the particle dynamics are 
given in the center-of-momentum (CM) frame where 
P = by the invariant mass operator, 



Mi = M + V. 



(2.15) 



The quantity M, introduced in Eq. (2.2), is the free in- 
variant mass in the Schrodinger picture and V is the el- 
ementary hadron interaction potential. 

The probability amplitude for observing an iV-body 
state |/3Q) with total momentum Q, given an initial N- 
body state \ocP) with total momentum P, is given by the 
^-matrix element (/3Q\S(E)\aP). The T-matrix T(E, P) 
is defin ed by the Lippmann-Schwinger equation (LSE) 
of Eq. (2.13) and determines the on-shell S'-matrix ele- 
ments, 



(f3Q\S(E)\aP) = (PQ\aP) 

-2m 5(8(Mp,Q) - E(M a ,P))(f3Q\T(E,P)\aP) 



(2.16) 



where 8(M a ,P) = y/Ml + P 2 . 

The potential V and the T-matrix describe all interac- 
tions between the various channels, including channels 
with differing numbers of particles. In general, they 
do not cons erve particle number. Therefore, the LSE 
of Eq. (2.13) represents an countably infinite system of 
coupled-channel equations which couple states of differ- 
ent numbers of particles. 

This infinite system of coupled equations may be sim- 
plified by truncating the Hilbert space to include only a 
finite number of states that are expected to contribute 
substantially to a given reaction. For the purposes of 
this study, the Hilbert space is restricted to contain a 
finite number of one-, two-, and three-particle states. 
Furthermore, the particles are assumed to be of finite 
spatial extension, thereby providing an ultraviolet reg- 
ularization to the theory. With these restrictions, the 
LSE of Eq. (2.13) reduces to a closed system of integral 
equations which may be solved exactly. Of course, one 



drawback of such a truncation is that some symmetries, 
such as crossing symmetry, which require the inclusion of 
many-particle states may be lost. The addition of states 
with a higher number of particles, such as four-particle 
states, can in principle be included in a straightforward 
manner but the resulting set of equations would be far 
more complicated than that studied here. 

The main objective of this work is to develop a frame- 
work for handling up to three-body channels in a fully 
unitary fashion, by including effects beyond their contri- 
bution to the real part of the effective, two-body poten- 
tials. The intended application is the description of soft 
final-state interactions in hadron production processes. 
Such processes are distinguished by their strong cou- 
plings and low momentum transfers. For this reason, 
composite hadrons (mesons and/or baryons) are chosen 
as the fundamental degrees of freedom rather than quarks 
and gluons. 

The truncation of the Hilbert space to contain only 
one-, two- and three- hadron states may be sufficient 
since, in many applications, states with higher numbers 
of hadrons contribute very little to two-hadron elastic 
scattering amplitudes. This suppression arises because 
many-hadron intermediate states typically have a large 
invariant mass, which appears in the denominator of the 
Green function G, tending to weaken its contribution. In- 
terestingly, this suppression of higher-order Hilbert space 
states is also observed in some quantum field theoreti- 
cal frameworks. In a study of the pion-loop contribtion 
to the p-meson self-energy and charge radius, based on 
a phenomenological application of the Dyson-Schwinger 
equations of QCD the covariant, quantum field the- 
oretic expression for the p-meson self-energy was sepa- 
rated into the various time orderings and their relative 
importance calculated. The time orderings include con- 
tributions arising from tttt and ppmr intermediate states, 
as well as others. In this calculation, it was shown that 
terms associated with the two-pion intermediate state 
contributed more than 95% of the total, while the four- 
hadron states contributed less than 5%. Thus, one ex- 
pects that a truncation scheme which neglects states with 
four or more hadrons should provide a reasonable descrip- 
tion of the residual strong interactions between mesons 
and baryons. 

The matrix elements of the potential V describe the 
couplings between hadrons that arise from the underly- 
ing QCD dynamics of quarks and gluons. Color confine- 
ment requires that all physical particle thresholds are as- 
sociated with the colorless hadron states. It follows that 
the matrix elements of the potential V are real. In this 
framework, all of the analytic structure of the T-matrix 
necessarily arises from the color-singlet hadron poles and 
branch cuts which result from the LSE of Eq. ( |2.13| ) . 

Once the Hilbert space has been truncate d to i nclude 
only one-, two- and three-particle states, Eq. ( 2.1 3| ) is ex- 
panded and rewritten in a simpler form by labeling each 
of the Hilbert space operators with subscripts indicating 
the numbers of particles they act on. The potential V is 



4 



V = 



> 




+ X 



o 



+ 







FIG. 1. Sch emat ic diagram of the interaction potential ma- 
trix V of Eq. (fll|). 



of the form 



V 



Vn V 12 V 13 

V21 V22 V23 
V31 V32 V33 



(2.17) 



The part of the potential associated with the coupling of 
a one-particle state to a two-particle final state is denoted 
V21. The resulting system of integral equations can be 
solved formally in a straightforward manner. 

It is important to note that each matrix element of 
the potential V in Eq. (2.17) is itself a matrix, since it 
may contain interactions between any number of different 
particle channels. That is, matrix elements of the form 
V21 describe the couplings of any one-particle state with 
any two-particle state. The number of one-, two-, and 
three-particle states one wishes to include depends on the 
specific application. For the application to 7T7t scattering 
considered in Sec. 
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a further simplification is made 
by assuming an absence of fundamental interactions in 
V connecting one-particle states to three-particle states, 
and three-particle states to three-particle states. Then, 
the potential V takes a simpler form, 



V n Via 

v= I v 21 v 22 v 23 
v 32 



(2.18) 



and is shown schematically in Fig. [|. 

The neglected terms Via = V 3 \ are associated with 
energy-independent transitions between one-particle and 
three-particle states. When such terms are neglected the 
only way in which a one-body state can decay into a 
three-body state is through a multiple-step process in- 
volving a two-body intermediate state. 

In setting the term V33 = 0, several possible elementary 
interactions have been neglected. First, V33 describes 
direct energy-independent couplings between two three- 
body states, as well as interactions in which two of the 
particles interact while the third particle is a spectator. 
Such terms may be important. One might argue that it 
is inconsistent to include direct two-body interactions in 
V22 , but neglect the analogous two-body (plus spectator) 
interactions in V33. Nonetheless, in this work such terms 
are ignored. The significance and role of these interac- 
tions will be addressed in future studies. 

In the truncated Hilbcrt space, the free Green function 
is a diagonal matrix 



II = 
K\ = X 
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FIG. 2. Schematic diagram of the two-particle se lf-ene rgy 
E, and the two-particle kernel K, as defined in Eq. (2.21). 



G = 



and the T-matrix is 



Gi 
G 2 
G 3 



(2.19) 



T = 



Til -M-2 Ti 3 

T21 T22 T23 
T31 T32 T33 



(2.20) 



In t he C M frame, each submatrix Gi, G2 or G3 in 
Eq. ( |2.19|) is itself diagonal since our hadron states form 
a complete, orthogonal set of eigenstates of the free in- 
variant mass operator M.. 

Up on in ser tion of the se fo rms for V, G a nd T from 
Eqs. ( pT8| ), ( plsj ) and ( fO(i| ) into the LSE fl2T3| ), one 
may formally solve this system of integral equations. 
It is convenient to consider the combination of terms 
V 2 2 + V23G3V32, which appears frequently in our formal- 
ism. These terms play an important role and so are col- 
lected and rewritten as the sum of S and K, 



Z + K =V22+V2 3 G 3 V 32 . 



(2.21) 



These are referred to as the two-particle self-energy S, 
and the two-particle kernel K . These terms are defined 
such that the matrix elements of the two-particle self- 
energy X contain only terms proportional to a <5-function 
in the relative momentum of the two-particle state. Con- 
sequently, matrix elements of the two-particle kernel K 
contain all contributions that are not proportional to a 
5-function in the relative momentum. 

The two-particle self-energy and kernel are depicted 
schematically in Fig. [| In the following, it will be- 
come apparent that S and K are the central elements 
of the framework, from which all other quantities are ob- 
tained. In fact, all effects due to three-particle interme- 
diate states can be traced back to these two amplitudes. 

One defines the dressed one- and two-particle Green 
functions in the usual manner as 



G 1 = (Gr 1 -n)- 1 , 

G 2 = (G 2 - 1 -S)- 1 . 



(2.22) 
(2.23) 



They are defined in terms of the two-body self-energy S 
and the one-body self-energy II, where 



n = v 11 + v 12 G 2 V2i, 
V21 = V21 + KG2V21, 

= (l - KG 2 )- 1 V 2 i. 



(2.24) 



(2.25) 
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FIG. 3. Schematic diagram of the integral equations for 
the one-body self-energy II and the dress ed ver tex V ix ■ These 
diagrams depict the expressions in (2.24) and (2.25). 
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FIG. 4. Schematic diagram of the scattering amplitudes 
t' 1 ^ and which enter into the two-particle scattering ma- 
trix eleme nts. Thes e diagrams depict the expressions in 
Eqs. (f2~27|) and (^28h. 



These quantities are shown schematically in Fig. ||. 

The solution for the two-particle scattering T-matrix 
can be written as 

T 22 = Gz 1 G 2 (t (1) + t^)G 2 G^ + G2 X G 2 E , (2.26) 



where 



t (2) 



— V 2 \G\V\ 2 , 
= K + KG 2 t (2 
= (I-KG2)- 1 



K. 



(2.27) 



(2.28) 



These two scattering amplitudes are shown schemati- 
cally in Fig. [|. Briefly, the contributions to the two- 
body scattering amplitude T 22 that proceed through one- 
body channels are denoted t^ x \ while contributions that 
don't proceed through one-body channels are denoted 
both iW and contain the effects of the two- and 
three-body singularities, but only contains one-body 
singularities. 

The matrix elements for the dressed Green function G 2 



defined by Eq. (2.23) are given by 



G 0a (p,E) = S afj (E-M ai2 (p)+ie)-Y lal3 (p,E) 



(2.29) 



One can collect the terms from Eq. (2.21) contributing 
to the two-body self-energy Y,p a , and organize them into 
the following sum, 



£ = <5S + V&GsVX' + VX>G 3 V, 



(i) 



r(2) 



(2.30) 

Here, the superscript (i) refers to the diagram in which 
the i th particle in the two-body state emits and subse- 
quently re-absorbs the particle 73. The term is iden- 
tified with the part of the potential V 22 that is propor- 
tional to a Dirac (^-function in the relative two- body mo- 
mentum. (All other terms that appear in Eq. ( 2.2l| ) but 



which do not appear in £ in Eq. (2.30), are part of the 
two-body kernel K .) 

U pon i nserting a complete set of three-body states into 
Eq. ( 2.3C ) and evaluating the resulting expressions in the 
overall CM frame with P = 0, one obtains 



/>oo 

„„„ Jo 



a 7i3( fc 13,P) 



x ^1713(^13) Gy(p,ki3,E) Ky 13Ql (fci3), (2.31) 



where 



G3V32 . A similar expression is ob- 
(2) 

tained for £^(p, E). The momentum integration is over 
the relative momentum k\ 3 between the first and third 
particles of the three-body intermediate state, J is the to- 
tal angular momentum of the system, the sum is over all 
three-body states 7, and the three-body Green function 
is 



G 7 (p, ki 3 ,E) 



1 



E-M 1123 (p, k 13 ) 



(2.32) 



For brevity the ubiquitous two-body phase space factor, 

P7i 3 ( fc 13) 



,(fcl3,P) 



k 2 



(2tt)3 2£(M^ 13 {k 13 ), P y 
and two-body Jacobian 

■M 7l3 (fci 3 ) 



P713 



13; 



2£ (m 7l , kis)£(m 73 , k 13 ) ' 



(2.33) 



(2.34) 



are introduced. The expression in Eq. (2.31) for the two- 
body self-energy Ej^p, E) is depicted in Fig. ||. The 
two-body self-energy Yip a (p, E) is then the sum, 

Zp a (p, E) = SE 0a (p) + Eg(p, E) + Sg(p, £0, (2-35) 
where the counter term is chosen to be 



ST, a (p) 



^(p,E) + ^(p,E)) . (2.36) 

E=M ai2 (p) 



This is necessary and sufficient to ensure unitarity and 
that the stable two-body system a± 2 is observed asymp- 
totically with the invariant mass M. ai2 {p). Evaluation 
of the two-body self-energy £ Q (p, E) requires calculating 
the i magin ary part and a principal part of the integral in 
Eq. (2.31). For energies E above a three-body thresh- 
old, these integrals encounter poles in the three-body 
Green function G 7 (p, ka.3, E) for values of the relative 
momentum ki 3 — k^, where fco 7 satisfies the relation 
M 1123 (p, k 0l ) = E. 

From Eq. ( j2.24| ), one obtains an expression for the one- 
body self-energy in the CM frame, 

1 f°° 
U a (E)=SU a + - V / dk a 7l2 (fc,0) 

Zm ^ 712 Jo 
xV aill2 (k) G 7 (k,E) Vy 12ai (k,E) . (2.37) 



G 




FIG. 5. Diagram depicting one of the contri butio ns to the 



two-body self-energy E^2(p, E) given in Eq. (2.31). Shown 



here, particle ai decays into particles 71 and 73, and these 
subsequently recombine to form particle j3\. The relative 
momentum between the intermediate particles 71 and 73 is 
fci3, the relative momentum of the incoming state |apip2) 
is p, and relative momentum of the outgoing state lyflqiqz) 
is q. Solid circ les denote matrix elements of the potentials 
from Eq. ( [2,18| ), V32 and V23, evaluated between two- and 
three-body states. 



In Eq. ( 2.37 ), m ai is the mass of the one-body state 
|aP), V aill2 (k) is the vertex function of the potential V12, 
Vy 12ai (k, E) is the vertex function for the dressed vertex 
V21, G y (k,E) is the dressed two- body Green function, 
and a 7l2 (&;,0) is a factor from Eq. ( 2.33 ) associated with 
the phase space of the two-body system 712. The one- 
body mass counter term 5H a is fixed by demanding that 
the elements of the one-body self-energies be identically 
zero when the driving energy E = £(m ai , P). In the CM 
frame, the mass renormalization condition is 



Re(n a (E)) 



E^m a 



= 



(2.38) 



In this framework, the finite size of the hadrons involved 
results in vertex form factors, such as Va l7l2 (k) , which fall 
off sufficiently rapidly with k to ensure the convergence 
of all integrals. Therefore, the counter terms oTI and <5E 
are both finite. 

Having obtained expressions for the dressed one- and 
two-boby Green functions, one next considers the two- 
body scattering amplitudes and t^ 2 \ These scatter- 
ing amplitudes depend on the one- and two-body Green 
functions, as well as the two-body kernel K . The two- 
body kernel is comprised of three contributions 



where 



K =K^+K^+K^\ 



V£pG 3 V 32 



(2) 



(2.39) 

(2.40) 
(2.41) 



and K^ 4pt ^ is the part of the potential V22 that is not pro- 
portional to a Dirac <5-function in the relative momentum. 
This latter term is the direct four-point coupling of four 
mesons, and it is depicted as four meson lines converging 
on a single point in Fig. [j] In Eq. ( [2.39 ), the parentheti- 
cal superscripts on the first two quantities refer to which 




FIG. 6. Definitions of the relative momenta and particle 



labels for the kernel K^(q,p) given in Eq. (2.44). 



of the two particles in the incoming state emits the ex- 
changed particle, the superscript on the third quantity 
refers to the direct four-meson interaction. In the first 
two terms, the exchanged particle is subsequently ab- 
sorbed by the other particle in the outgoing state. 

Explicit ex press ions f or the most general two-body ker- 
nels of Eqs. ( 2.40 ) and ( 2.41 ) in the spherical wave basis 
are complicated and not particularly enlightening. How- 
ever, in the model application to tttt scattering consid- 
ered in Sec. Ill, the resulting kernel is relatively sim- 
ple. The only matrix elements of V23 that are of interest 
in this application are those associated with the transi- 
tions of the forms irn — > intp, tttt — * tttt/q, tttt — > KKp 
and tttt — + KKfo, tttt — > ttKK*. In each of these 
hadron states, at least two of the three particles are spin- 
mesons. For these interactions, the plane- wave matrix 
elements of the potentials are of the form 

( 7l 23Kkki3|^ 3 ( 2 1) |ai2Pp) = (2tt) 3 <5(K - P)<5 72 , Q2 

x2E (m a2 , -p) 
x(7 13 kk 13 |^2i|aip), (2.42) 

where the matrix elements of V12 are of the form 

(-1)^3 ~ A 



(7i3kki 3 |V2i|aip> = (27r) d J(k - p) 



73 



1 



X^5,A(-ki3.p) Ki* r s ; si _ A (k 13 ) v 113ai (k 13 ). (2.43) 

A 

The vertex for V32 appears in Figs. [| and [|as the right- 
most interaction vetrtex. In the partial-wave basis, the 
kernel is shown in Fig. ^ and is given by 



K^ J ( q ,p)=2n 



7123 



1 





1 47r 2c(to 73 ,p 



q) 



^2723 ( fc 23) S s ™ (q,p) V 113ai {kiz) 

E - ^7123 (q. -p. p - q) + ^ 



(2.44) 



where 



S s -= u (q,p) 
^3=l( q , p ) 



<4 



fc 73Al fc 7 3„ 



•M 723 (fc23)Al7l3 ( fc 13) 



(2.45) 
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for scalar exchanges (s 73 = 0) and vector exchanges 
(s 73 = 1), respectively, and Pj(x) are the usual Legen- 
dre polynomials in x — p-q/pq. The three four- momenta 



appearing in Eq. (2.44) are 



<4 = (£K3!,q),q) , 

(£(m 73 ,p-q),p-q) 



^73 



(2.46) 



Expressions corresponding to the matrix elements of 

Kpail'P) can ^ e obtained in a similar manner. 

Once specific forms of the model vertex form factors 
Kn^i (/C1 3) an d V^ 2723 (fc23) are provided and substituted 
into Eq. (2.44), the matrix elements of the kernel K are 
computed numerically. One can proceed to solve the in- 
tegral equation in Eq. ( 2.28| ) for the scattering ampli- 
tude . In the CM frame, the integral equation for the 
partial- wave scattering amplitude (q>p) has the form 

poo 

t?l J (Q,P) = K^ a (q,p) + Y^ / dka yi2 (k,0) 

712 Jo 

xK^(q,k) G 7 (k,E) t^ J (k,p), (2.47) 

where a 7l2 (fc,0) is the usual two-body phase space fac- 
tor. Obtaining the solution of this integral equation is 
complicated by the presence of the two-body pole in the 
two-body Green function Gj(k), and possibly the appear- 
ance of three-body unitarity cuts in both G 7 (k) and the 
two-body kernel Kp a (q,p). The method used to solve 
this integral equation is adapted from Ref. p0[. It in- 
volves obtaining a two-body Moller operator fK^ , whose 
J th partial-wave matrix element satisfies, 



712 



tpl J (q,p) = T,L dk «7i 2 (M) n$ J (q,k) K ja (k, P ). 

(2.48) 



After a solution for t^ J (q,p) is obtained, one proceeds 
to obtain an explicit expression for the two- body scatter- 
ing amplitude which is given by Eq. ( 2.27 ). Since 
the solution of the intermediate dressed one-body Green 
function G\ is obtained from the one-body self-energy 
from Eq. (2.37), all that remains is to determine the form 
of the dressed vertices V\2 and V21. The integral equa- 
tion for the dressed vertex is o btain ed for the transpose 
of the dressed vertex from Eq. (2.25), 



1*00 

Kaisai (<?) = ^120! (?) + / * <h™ 
712 Jo 



(fc,0) 



rJ = S 



x K" 01 Q1 (q, k) G 7 (fc, E) V 112ai (k) (2.49) 

The similarity between t his in tegral equation and the in- 
tegral equation of Eq. ( 2.47 ) with J = s ai which de- 
termines is clear. It follows that the solution to 



V /3i 2 <*i (?) = / dk a 7l2 ( k , 0) ^al J (<?: k ) V ll2ai (k). 

7!2 J ° 



(2.50) 

Finally, the two-body scattering amplitude is 



7i7i 



2^/TO 7 jTO 7l 



The complete expression for the two-body scattering am- 
plitude T22 is obtained b y add ing this expression for 
to according to Eq. ( 2.26 ). 

In the previous sections, it was demonstrated that the 
explicit solution to the scattering problem involving one- 
, two-, and three-body states can be obtained by per- 
forming several integrations and one matrix inversion. 

The matrix inversion is necessary to obtain the two-body 

(2) j 

Moller amplitude fi^ (q,p). 



III. APPLICATION TO tttt SCATTERING 

In this section, the framework is applied to tttt scatter- 
ing. Simple model forms of the elementary vertex form 
factors Vf3 12ai (q) are introduced, and solutions for the 
self-energies Hp a (E) and T,p a (p,E) and scattering am- 

(2.) 7 ~ 

8a (9>-P)i and Vp 12ai (q) are obtained 



plitudes 4^(9,2/), t 
numerically. Several interesting aspects of the obtained 
solutions are discussed. It shoul d be emphasized that 
the model introduced in Sec. Ill A is preliminary and the 
manner in which the model parameters are fit to the data 
may be overly simplistic, as it focuses on reproducing 
only a few observables and therefore does not represent 
an exhaustive or complete study of the dynamics of tttt 
scattering. The motivation is to provide a demonstration 
of the framework and exhibit the features of the model, 
and its ability to describe the scattering of a system of 
strongly-coupled particles with emphasis on the multi- 
particle channel aspect. More complete studies of meson 
scattering within the present framework will be the sub- 
ject of future articles. 

In Sec. Ill A, the dynamical assumptions are discussed 



Eq. (2.49) is just 



along with the model parameters. A detailed list of the 
states included in the Hilbert space is provided. The 
model parameters are determined using a simple method 
to fit experimental data for the 7T7t isoscalar-scalar phase 
shift and p-meson decay width ; usin g the S'-wave phase 
shifts from Ref. In Sec. [II B the resulting phase 

shifts, inelasticities and cross sections are provided and 
compared to the data, and some aspects of the /o(980) 
scalar meson are discussed in terms of a K K bound state. 



A. Dynamical model for nn scattering 

The model is intended to describe the scattering in a 
range of center-of- momentum (CM) energies from thresh- 
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old (E = 2m^ w 280 MeV) to about E = 1400 MeV. 
Above 1400 MeV, it is important to include in more de- 
tail the effects of the three scalar mesons observed in this 
region. For this preliminary study, however, it is possi- 
ble to avoid making strong assumptions concerning these 
scalar mesons, hence the model will not be accurate in 
this energy region. In the following, only the isoscalar- 
scalar I = 0, J = 0) and isovector- vector (I = 1, J = 1) 
channels are considered. The motivation is to explore 
some of the interesting physical aspects of the present 
framework and to estimate the importance of including 
three-body states in such a model of hadron scattering. 
The assumptions of the dynamical model are summarized 
below. 

Two-body states: For the channels and energies ex- 
plored herein, it is assumed that tttt scattering is primar- 
ily determined by the dynamics arising from the coupling 
of the 7T7T and KK two-body channels. Hence, |tttt) and 
\KK) are the only two-body channels included in the 
Hilbcrt space. 

One-body states: It is assumed that the coupling of 
the KK system is strong enough to result in the appear- 
ance of a narrow resonance in the scalar-isoscalar chan- 
nel at E fn 980 MeV. This state is identified with the 
J PC = ++ /o(980) meson. Since this scalar meson is 
presumed to arise from final-state interactions as a quasi- 
bound KK state, it is not part of the free Hilbcrt space, 
and there is no bare mass associated with it. Rather, it 
appears as a pole in the analytically-continued T-matrix. 
Furthermore, in the limit that the two-body tttt and KK 
channels decouple, this pole moves to the real-energy 
axis below the two-kaon threshold; that is, it becomes 
a KK bound state in this limit. The identification of 
the /o(980) meson as a KK molecule is controversial. 
Although it appears as a molecular state in this model, 
the "true" nature of the /o(980) meson remains an open 
question. 

In contrast to the /o(980) meson, it is assumed that at 
least one of the scalar resonances observed in the mass 
region between 1300 and 1700 MeV will be a QCD bound 
state; that is, a state which arises as a bound state whose 
constituents are quarks, antiquarks and gluons. Such 
states do not arise from the meson final-state interac- 
tions, they are not bound states of mesons, and hence 
must be included in the model as bare states with bare 
masses. 

Experiments reveal the presence of several reso- 
nances in the scalar-isoscalar channel between 1300 and 
1700 MeV. A complete study of the tttt scattering sys- 
tem in this energy range requires the inclusion of each 
of these resonances into the model. However, to simplify 
the present study, all of these resonances are modeled 
in terms of a single scalar resonance. The resonance is 
assumed to have a mass of 1350 MeV, which gives it 
a mass similar to the lightest of the resonances above 
kaon threshold, referred to as the J PC = ++ /o(1370) 
meson. One ramification of choosing a single scalar res- 
onance to model the effect of all observed resonances in 



TABLE I. Masses and widths of mesons. Masses that are 
underlined have been fixed to reproduce the accepted values. 
All other values are obtained from the model calculation. In 
the present study, the width of the K* meson was not calcu- 
lated. All values are given in units of MeV. 

it K p K* /o(1350) /q(980) 
mass 140 500 770 890 1350 996 
width 150 805 46 



TABLE II. Isospin coupling constants for kernel Kp a {q,p) 
for isoscalar, S-wave (7 = 0, J PC = ++ ) and isovector, 
P-wave (I = 1, J PC = 1 ) scattering. 



Channel 


Exchange 


1 = 0, J = 


1=1, J = 1 


7T7T <-> 7T7T 


tttt/o 


1 


1 


KK <-> KK 


KKf 


1 


1 


Tin <-> 7T7T 


imp 


-1 


-1 


KK <-> KK 


KKp 


-1 


-1 


Tin <-» KK 


■kKK* 


-V2 


-1 



the 1300-1700 MeV region, is in the width of this reso- 
nance. In order to fit the model parameters to the tttt 
phase shifts requires a single effective resonance with a 
very large width. It is found that the model resonance 
has a decay width of 805 MeV, which is approximately 
the sum of the widths of the three observed resonances 
in this region. 

Three-body states: Only three-body states are included 
that can couple to the 7T7t or KK states through the ab- 
sorption or emission of the isovector J PC = 1 p(770), 
isodoublet J p = 1" if* (892), and J PC = 0++ f mesons. 
Thus, the three-body states included in this study are 
\rnrp), Itttt/o), \KKp), \KKf ), \ttKK*) and \nKK*}. 

To summarize, the hadronic states included in this 
model of 7T7t scattering are 

l/o>, \P), 
M, \KK), 

\wnp), \KKp), \ttKK*), Itttt/o), \KKf ), (3.1) 

where the /o meson refers to the /o(1350) meson. (The 
/o(980) meson is expected to appear in the model as a 
KK resonance.) The values of the bare masses of these 
particles are provided in Table | and are underlined to 
indicate that they are input parameters. As discussed 
m Sec. ||, one- and two-body counter terms are included 
in the elementary interaction potentials V\\ and V22, re- 
spectively, such that the bare masses given in Table || 
coincide with the dressed masses of the mesons. 

Model vertex form factors: The vertices in the model 
are assumed to be finite-sized and hence require the ap- 
propriate form factors for the relative three-momentum 
q. They are given by the universal form: 



Vp 12ai (q) = a PlP2ai J e-" 2 /A| lfcQ1 , (3 . 2 ) 

V "omax ~r 1 
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where s max = max{s J 3 1 ,S/3„s ai } is the largest spin of 
the particles involved. In the present study s = for 
vertices involving the fo meson, and s = 1 for ver- 
tices involving the p or K* mesons. The vertex cou- 
pling constants a / 3 1/ g 2Ql and form factor momentum scales 
A^j^qj are chosen to provide a good fit to the data 
for the isoscalar-scalar rnr phase shift 5^ and the p- 
meson decay width T p ^,^ = 150 MeV. The parameter 
search was limited in a number of ways. First, the var- 
ious meson-exchange form factor scales A^ 1 ^ 2ai in Ta- 
ble [II were all constrained to be the same value and 



less than 1 GeV. The direct four-meson couplings in Ta- 
ble IV were chosen to be one of two scales, the first 



was taken to be 125 MeV larger, and the second to be 
125 MeV smaller than the meson-exchange scales in Ta- 
ble III. In the following, it is shown that the vector- 



exchange interactions contribute little to the observables 
considered. To reduce the number of parameters, the 
strength of the vector-exchange vertices were taken to be 
identically equal a^p = ciKKp = a vKK* ■ The isospin 
factors that arise in a calculation of the meson-exchange 
kernels Kl a (q,p), such as in Eq. ( jUg ), are given in Ta- 
ble n 

Direct interactions: In addition to meson exchanges, 
it is important to include real-valued potentials that 
directly couple two pseudoscalar mesons to two pseu- 
dosc alar mesons, as a part of the if( 4pt ) kernel in 
Eq. ( 2.39 ). Such interaction potentials could arise from 
the direct coupling of four mesons to a virtual-quark loop. 
Here, two direct four-point interactions are considered. 
The first is intended as a way to mimic some of the ef- 
fects of dynamical chiral symmetry breaking. This inter- 
action is taken to be of the form given by the elementary 
potential V22 and is referred to as the direct 4tt (or 4K) 
interaction. In a partial-wave basis, the form of the four- 
pion interaction is given by 

Kt£(q,p) = 167T (qpY a% e"^ e^L (3.3) 

and the four-kaon term K^ 4K ^ J (q,p) = 0. The second 
four-point interaction is a short-ranged attraction mod- 
eled as a i-channel exchange of a heavy scalar-isoscalar 
meson. Its form is given by the scalar-exchange kernel of 
Eq. (|2.44[) and the two-body self-energy S of Eq. ([2.3l[) . 



For simplicity it is treated exactly as if two additional 
three-body states, 



IttttV), \KKX) 



(3.4) 



with mx = 1500 MeV, were added to the Hilbert space. 
Clearly, modelled in this manner, for energies E > mx + 
2m OT = 1780 MeV, the state \mrX) can go on-energy- 
shell. However, the calculations described herein are for 
energies less than 1400 MeV, so that \ttttX) is never on- 
energy-shell. 

Again, the objective is to study the framework devel- 
oped in this paper. That is, it is interesting to assess 
the importance and possibility of including the dynamics 



TABLE III. Coupling strengths and momentum scales for 
the vertex form factors. The values of the momentum scales 
are given in MeV. 



7T7r/o KKfo 7T7T/9 KKp K-kK* 



A/3i/3 2 oi 



12.4 
875 



5.08 
875 



20.0 
875 



20.0 
875 



20.0 
875 



tyttX 


KKX 


4n 


4K 


19.8 


11.0 


17.2 


0.0 


1000 


1000 


750 


750 



TABLE IV. Coupling strengths and momentum scales for 
the four-point meson interactions. The values of the momen- 
tum scales are given in MeV. 



a 
A 



of three-body intermediate states to the study of meson 
scattering and final state interactions, rather than to test 
a particular interaction model for ttt: scattering. 

Once the forms of the model vertices are fixed, the 
only observables used in the fit are the scalar-isoscalar 
7T7T phase shift S^, the existence of the K K resonance 
(referred to as the / (980)), and the decay width of the 
p(770) vector meson. All measurable reaction channels 
are not used in the fitting procedure, since this paper 
represents more a proof of principle of the framework 
rather than a complete phenomenological analysis of 7T7t 
scattering. The resulting values of the coupling constants 



are provided in Tables [II and IV 



B. Phase shifts and inelasticities 

Below all three-body thres holds, the non-trivial part 
of the S'-matrix in Eq. ( 2.16 ) can be written as a 2 x 2 
unitary matix Sp a (E) with a,0 denoting the only two 
open channels 7T7t and KK. The S-matrix in the J th 
partial wave can be parametrized in terms of two phase 
shifts Smr and S K g, and one inelasticity T)^, 



2 r)a(E) 
Va(E) = \S J aa (E)\, 



(3.5) 
(3.6) 



for a — 7T7T, KK. Below all three-body thresholds there 
are only two stable channels, irir and KK. Hence, there is 
only one inelasticity parameter r]^ = r\ K ^ . For energies 
E above the lowest stable three-body threshold, one must 
augment the S'-matrix by including all stable three-body 
states. Consequently, its parametrization requires more 
than two phase shifts and one i nelas ticity. Nonetheless, 
one may still use Eqs. (3.5) and (3.6) to define the phase 
shifts 5 a (E), and inelasticities rj a (E) for the two channels 
a = 7T7T and KK. Of course, above the threshold of a 
stable three-body state 7/ ff7r 7^ Vkk- 

The 7T7T phase shift, as defined by Eq. ( |3.5| ) and ob- 
tained from our model, is shown as a solid curve in Fig. [7]. 
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0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 

E (GeV) 

FIG. 7. The tt-k scattering phase S^n in the scalar, isoscalar 
channel as a function of the CM driving energy E. The data 
are from Refs. [^2j (open squares), |^| (closed circles), |2j|] 
(up triangles), and p5| (down triangles). 



The model provides an excellent description of the 7T7r 
phase shift data depicted in Fig. fj] [p^-p5|, and the in- 
elasticity ?7 W7r , as shown in Fig. ||. The overall trend of 
the pion scattering phase shift 5 VV is positive and in- 
creases slowly with energy E. This is indicative of a weak 
and attractive effective 7T7t scattering potential. At the 
kaon threshold, a rapid phase motion is apparent. It re- 
sults from the presence of a narrow, /o(980) scalar meson. 
Above the two-kaon threshold, the phase shift continues 
to increase slowly, at a rate similar to the increase in the 
phase shift below the threshold. 

Below all other thresholds, the 7T7t channel is the only 
open channel, and unitarity requires that the inelasticity 
ifrn = 1 here. This is clearly observed in Fig. ^, where 
the calculated inelasticity rj^ has a value consistent with 
unity below the threshold of the KK channel at 1 GeV. 
Had four-body states been included into this framework, 
one might have expected to see a decrease in the inelas- 
ticity ^tj-tt due to the opening of the four-pion state, which 
has a threshold of E = Am n sa 0.560 GeV. However, the 
data in Fig. ||| |22| seem to suggest that the contribution 
of the four-pion state to 7T7t scattering is negligible. This 
can be seen by noting the lack of any systematic deviation 
from rjim = 1 for the range of energies im,^ < E < 2m#-. 

It is clear from both Figs. [?] and || that the KK channel 
has a significant effect on 7T7t scattering. At the two-kaon 
threshold at E w 1.0 GeV, one observes a rapid increase 
in the 7T7t phase shift , and a sharp fall off of the 7T7t 
inelasticity rj^ to its minimum value rj vn w 0.31. This 
phase motion is indicative of crossing the thresholds of 
the two-body KK state and the one-body /o(980) bound 
state. The rapid increase observed in the phase shift 
i5 W7r (E) is due to the weak coupling between the 7T7t and 
KK channels. In the model, when the mixing of these 
two channels is further weakened, the rate of change of 
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FIG. 8. The hit scattering inelasticity r\^^ in the scalar, 
isoscalar channel as a function of the CM driving energy E. 
The data are from Ref. p3|. 



the phase motion tends to increase, until finally, in the 
limit that the coupling between the two channels goes to 
zero, the phase motion becomes a step-function of mag- 
nitude 180 degrees. Such a phase motion is completely 
unobservable, and could therefore be ignored altogether 
(although it is relevant to Levinson's theorem, which re- 
lates overall changes in the phase shifts from threshold 
to infinite energy to the number of bound states in a 
system). 

The importance of the coupling between the two-body 
channels 7T7t and KK can be estimated quantitatively by 
recalculating the pion phase shift 5^ after removing the 
kaon state \KK) from the Hilbert space. The result is 
shown as the dotted curve in Fig. |^. However, from the 
above argument, perhaps a better indication of the im- 
portance of the KK channel is obtained by letting the 
couplings that lead to a mixing between the KK and 7T7t 
channels go smoothly to zero. (In practice, this is done by 
not allowing the one-body scalar /o(1350) to have a bare 
coupling to the KK state, and setting a,K*nK — 0. This 
has a minimal effect on the dynamics but prevents mixing 
the |7T7r) and \KK) states.) In this limiting case, the two 
states \KK) and |/o(980)) are coupled to the |7T7t) state, 
but the contributions they make to 7T7t scattering go to 
zero. The result is that for energies E > 1000 MeV, the 
7T7T scattering amplitude will cross the branch cuts asso- 
ciated with the opening of these two channels. But since 
these channels do not mix with the 7T7t channel in this 
limit, the 7T7T phase shift exhibits a step-like motion at 
E = 1.0 GeV. This is shown as the dashed curve in Fig. ||. 
The difference between the solid curve, which represents 
the full model calculation, and the dashed curve in Fig. ^ 
may be taken as being indicative of the significance of 
the KK channel on 7T7t scattering. One concludes that 
the mixing between the 7T7t and KK channels is signif- 
icant near and below the two-kaon threshold, where it 



11 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 

E (GeV) 

FIG. 9. The tt-k scattering phase S^n in the scalar, isoscalar 
channel versus the CM energy E. The full calculation (solid 
curve) is compared to the calculation (dashed and dotted 
curves) in the limit the KK and irn states become decou- 
pled. The data are from Ref. pi. 



can contribute more than half of the total phase shift 
d n7r . At energies above the two-kaon threshold, its im- 
portance quickly dimishes and vanishes altogether above 
E = 1150 MeV. 

The importance of the mixing between the itit and KK 
channels can also be observed in the pion inelasticity rj^ 
shown in Fig. ||. Just above the two kaon threshold, the 
inelasticity plummets to a minimum value rj^^ m 0.31. In 
the limit that the coupling to the kaon channel goes to 
zero, as described above, the inelasticity takes on a very 
different appearance. 

In Fig. 10, the two pion inelasticity is plotted above the 
KK threshold. Below the threshold its value is unity. 
The lowest multi-particle state to which two-pion flux 
can be lost is the three-particle channel \mrp). The pro- 
duction threshold of this state is 1050 MeV. It is clear 
that one observes a slow decrease in the inelasticity 
above 1050 MeV, due to the opening of the 7T7rp channel. 
The effect of this channel on the inelasticity is very small, 
with a minimum value 77^ w 0.994. 

The energy dependence of the S-wave pion scatter- 
ing inelasticity 77^ and phase shift S ww are conveniently 
plotted together in an Argand diagram. In Fig. the 
partial- wave scattering amplitude a^°(E) is plotted in 
the complex plane as a parametric function of the CM 
energy E. This figure clearly shows the rapid rise in the 
elastic irir phase shift just below 1 GeV, and the resulting 
dramatic loss of flux from the elastic channel as soon as 
the KK channel opens up. 

For spinless particles, the amplitude is given by 
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FIG. 10. The irn inelasticity 77,^ in the isoscalar-scalar 
channel with no kaons. 
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x fei>o,Po) + t£>i»(Po,Po)) , (3 



where po = \J E 2 /4 — m 2 ^ is the magnitude of the on- 
energy-shell three-momentum of the pions, 07T7T is the 
wave function renormalization of the two-pion state |7T7r), 



and t^ J (q,p) and t^ J (q,p) are the two-body scattering 



AWi 



amplitudes, obtained from Eqs. (2.47) and ( [2.51 ), respec- 
tively. The two-body scattering cross section a/3^ a (E) 
can be written in terms of the partial-wave scattering 
cross sections according to 



ap^ a (E) 



J=0 



(3.9) 



In terms of the scattering amplitudes, these partial-wave 
cross sections are given by 



4^{E) = 
2 J + 1 zpz a 
4tt Mtt 2 E 2 



t { ) a J (qo,Po)+tf a \qo,Po) , (3.10) 



where po and go are the on-energy-shell solutions to 
■M. a {po) = E and Mp(qo) = E, respectively. The re- 
sulting cross section for elastic, S- wave 7T7t scattering is 
shown as a solid curve in Fig. [l2]. It is finite at thresh- 
old, exhibits a maximum value of 43 millibarns at E « 
600 MeV, and a sharp decrease at the position of the 
/o(980) scalar meson resonance. This sudden drop oc- 
curs just below the KK threshold, as can be seen upon 
examination of the inset plot in Fig. [l2] which depicts a 
closeup of the KK threshold region. 

The dot-dashed curve is the resulting cross section 
cr^P^^ for the two-kaon production process 7T7t — > KK. 
This cross section is considerably smaller than that of the 
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FIG. 11. Argand diagram for nn scattering at all energies. 
Plotted is the energy dependence of the real and imaginary 
parts of the partial- wave scatte ring amplitu de aj^° (E) for 
the J = partial wave from Eq. (3.7) or Eq. (3.8). The curve 
is calculated for all energies E > 2m- K , and annotated with 
the corresponding energies E in GeV, from threshold to above 
E = 1.7 GeV. 



elastic rnr scattering cross section, reaching its maximum 
value of 4.6 millibarns just above the two-kaon threshold 
at E = 1 GeV. Its small size is a result of the weak cou- 
pling between the two-pion and two-kaon channels. As 
discussed above, a weak coupling of these channels is nec- 
essary to ensure a narrow /o(980) meson. If the mixing 
between the pion and kaon channels were stronger, the 
/o(980) meson would more easily decay into two pions, 
tending to increase its width significantly. 

The dashed curve in Fig. [l^ is the S- wave cross section 
for elastic KK scattering. This cross section is compara- 
bly huge, having its maximum at the two-kaon threshold 
energy. Its size can be compared to that of the two-pion 
elastic cross section, 



a J=o 

7T7T< — 7T7T 

T J=P _ 



(E = 2m„) w 18.3 mb, 
(E = 2m K ) « 734 mb. 



The very large KK cross section arises from the scalar 
/o(980) meson which lies just below the KK. The pres- 
ence of a bound state just below the two-body scattering 
threshold will generally tend to increase the size of the 
cross section dramatically. 



As was discussed earlier in Sec. Ill A, the model pa- 
rameters were chosen to provide a good fit to the S'-wave 
pion phase shift S„„ data from Ref. plj . In Fig. [7| the 
resulting phase shifts are also in excellent agreement with 
the data analyses of Ref. p^ , p4| ], and p5| ]. However, 
it is important to realize that direct comparison of these 
model results with these nir scattering data must be done 
with caution. Extraction of these data requires the use 
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FIG. 12. S'-wave cross section a J=0 (E) in millibarns for 
7T7T — > 7T7T (solid curve) , KK — > KK (dashed curve), and 
7T7T — ► KK (dot-dashed curve). The inset shows a detail of 
the region around the KK threshold at E = 1 GeV. 



of theoretical models, or theoretical assumptions, in or- 
der to fit the experimental observables. In the worst-case 
scenario, the resulting 7T7t phase shift data may be more 
representative of the extraction methods employed than 
of the actual irir scattering process. 

The extraction of the 7T7r phase shifts from experi- 
ment is a difficult and long-standing problem of hadron 
physics. At present, it is impossible to construct an ex- 
periment in which a pion beam is scattered from a pion 
target. Hence, other techniques are required to extract 
the 7T7T phase shift from experimental observables. One 
possibility is to use decays that produce two pions in the 
final state, and attempt to extract the 7T7t phase shifts 
from the final-state interactions. 

The procedure employed by Ref. |24|] is extract the 
phase shifts from the electroweak kaon decay K + — > 
ir + ir°e + v e . The results are shown as open squares in 
Figs. [)]and E3[ Extraction of the phase shift using decays 
with more than two particles in the final state requires 
some knowledge of transition form factors for the cou- 
pling of a kaon, two pions and the W boson to determine 
K + — » ir + ir~W + . A nice feature of employing the elec- 
troweak decay K + — > 7r + 7r°e + ^ e is that the two pions are 
the only strongly-interacting particles in the final state. 
Hence, one would expect that the 7T7t interactions would 
be the dominant contribution to the dressing of the final 
state. However, this approach is hampered by two ex- 
perimental difficulties. The first is the lack of statistics. 
This particular kaon decay represents only a small frac- 
tion (4 x 10~ 5 ) of the total K+ -meson decay width, which 
is already extremely small. The second is the fact that 
the energies and angles of the outgoing leptons provide a 
small lever arm with which to vary the CM energy E of 
the two final-state pions. 
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Another method is to extract the final state interac- 
tions of the two-pion production process irp — > mrn, 
employed by Ref. [§5l (down triangles), Ref. [^2| (open 
squares), and Ref. |23[ | (closed circles), shown in Fig. |j. 
These studies require some theoretical input in order to 
perform the extraction of the tttt scattering phase shift 
8-ktt and inelasticity 77^; hence, they are not direct mea- 
surements of the 7T7T scattering phase shift. 

Our model parameters were originally fit to the mr 
phase shifts obtained by an analysis |23|| of an exper- 
iment at CERN involving ir~p — > 7r + 7r~n at 17.2 GeV. 
Recently, this same data was re-examined by Kaminski et 
al. p2| , with weaker model assumptions than were used 
in Ref. [^3| . The work of Ref. |^2| provides an exhaus- 
tive and nearly complete study of the mr phase shift. In 
particular, there is no assumption that pion exchange is 
the dominant mechanism for the process ir~p — > ir + Tr~n. 
Consequently, this analysis seems to be more general 
than the others. A relative phase ambiguity in the anal- 
ysis of Ref. |E| provides four possible, distinct solutions 
for the 7T7T phase shifts and inelasticities. Two of these so- 
lutions seem to have an unphysical inelasticity tj vv below 
the two-kaon threshold and can be discarded. The other 
two solutions, denoted the "up-flat" and "down-flat" so- 
lutions, are very similar in appearance and neither can 
be dismissed on qualitative grounds. The phase shifts 
from the "down-flat" solution of Ref. 0] are shown as 
solid circles in Figs. || and |l3|. 

The general behavior of the mr phase shift 6^^ is pos- 
itive, which is indicative of an attractive mr scattering 
amplitude. An abrupt increase in the pion phase shift is 
evident at 1 GeV, which is due to the combined effects 
of the opening of the KK threshold and the crossing of 
the scalar /o(980) resonance. 

Apart from this feature, which in this particular model 
is a result of the delicate mixing between the mr and 
KK channels, the calculated pion phase shift <5 T7r ex- 
hibits a steady, gentle increase from the threshold at 
280 MeV to above 1400 MeV! In the model, this slowly 
increasing behavior arises from subtle cancellations be- 
tween the attractive potentials of the heavy scalar-meson 
exchanges and the repulsive, direct four-pion interaction 
from Eq. (3J3) which is intended to model the effect of 
dynamical chiral symmetry breaking. The different form 
factor scales involved in these interactions (see Tables [II 



and IV) are chosen to provide this slowly increasing, weak 
phase shift observed in the two-pion channel. Typically, 
scalar potentials by themselves provide a strong attrac- 
tion in the 5-wave mr channel, that leads to a rapidly 
rising phase shift just above the mr threshold, which then 
quickly falls away. This behavior is not seen in the mr 
phase shift 5 n7r . 

The importance and role of the scalar resonances in 
mr scattering can be appreciated by a close examination 
of Fig. [l3|. The effect of the heavy s-channel resonances 
which are collectively modeled by the single /o(1350) in 
the model is two fold. First, their presence leads to a 
strong attractive potential for energies below 1400 MeV. 



360 

315 

270 

225 

?180 
■a 

^135 
to 90 

45 




1 1 1 1 1 1 1 II 1 1 II 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 II 1 1 II 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

1 1 1 — 




"1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 f 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 

E (GeV) 

FIG. 13. Importance of the different model contributions 
to 7T7T scattering phase shift. The solid curve is the full calcula- 
tion. The dashed curve is obtained by removing the one-body 
state /o(1350) from the Hilbert space. The dot-dashed curve 
is obtained by completely removing the scalar /o(1350) from 
the theory, both in one-body and three-body states. The 
dotted curve is same as the dot-dashed curve, but all scalar 
couplings (the "X" -exchanges and contact terms) are also set 
to zero. This is the effect of just the vector-meson exchanges. 
Above threshold it is negative with a minimum of S^n ~ -2.0 
degrees at E m 1.2 GeV. The dot-dot-dashed curve is ob- 
tained by weakening the coupling of the one-body /o(1350) 
state to the irir state by 10 percent. The data points are from 
Fig. §. 



Second, they provide the most important contribution to 
off- diagonal matrix elements of the two-body scattering 
kernel K. That is, they provide the strongest source 
of mixing for the two-pion and two-kaon states in this 
model. 

Both of these effects tend to produce an attraction for 
the two pions. In particular, the strong attraction nec- 
essary to bind the kaons to form the /o(980) resonance 
results in a strong attraction in the two-pion channel as 
well. The amount of mixing between the two-pion and 
two-kaon states dictates the attraction felt by the pions. 
Hence, when the coupling between |7T7t) and |/o(1350)) 
states is artificially reduced by as little as 10%, the re- 
sult is significant, as shown by the dot-dot-dashed curve 
in Fig. [yj. When the one-body state |/o(1350)) is re- 
moved entirely from the Hilbert space, the result is the 
dashed curve. The resulting pion phase shift is negative 
and close to zero below the two-kaon threshold, and pos- 
itive above the threshold. The absence of the one-body 
s-channel state /o(1350) reduces the mixing between the 
two-pion and two-kaon states, which results in a nearly 
stable (and very narrow) kaon bound state Jo (980). In 
this case, one observes a Jo (980) bound state with a width 
that has been reduced from 46 MeV to 0.28 MeV! This is 
a result of the fact that the fo — > mr decay must proceed 
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through K* exchange in the kernel K, which provides 
only a weak mixing of the tttt and KK states. 

When all scalar mesons are removed from the the- 
ory entirely; that is, when the couplings that lead to 
the existence of one-body states |/o(1350)), and three- 
body states |/o7T7r) and \/qKK) are set to zero a^fa = 
a KRf = 0; the resulting phase shift is shown in Fig. [l3| 
by the dot-dashed curve. The slightly repulsive behavior 
is a result of the combined effect of the attractive scalar- 
X and repulsive four-point interactions associated with 
chiral symmetry given by the parameters in Table IV. 



To minimize the number of free parameters all of the 
couplings to vector mesons are chosen to be equal to each 
other a^-Kp = o-KKp = a irKK* ■ These coupling strengths 
were then determined by solving P-wave tttt scattering 
at the p meson mass E — m p — 770 MeV, and requir- 
ing that the p-meson width reproduced the experimental 
value r p = 150 MeV, as given in Table H. It is found 
that the resulting coupling strength leads to a vector- 
meson exchange interaction kernel K which provides a 
very weak repulsion for tttt scattering. This is illustrated 
by the dotted curve in Fig. |l3|, where all of the cou- 
plings except those involving the vector mesons (a p7r7r , 
a pKKi a K*-KK) are se t to zero. The resulting phase shift 
is negative (repulsive) and very small; its largest absolute 
value is about 2 degrees. Thus, in this model, p-meson 
exchange is negligible in the S wave. This model differs 
from the analysis of Ref . JlJ] , in which they report that 
the attractive potential (which leads to the binding of 
the KK into the Jo (980) resonance) is primarily due to 
p-meson exchange, which is strong and attractive in their 
model. In the framework, the exchange of a spin-1 meson 
in the kernel K results in a very weak and mildly repul- 
sive interaction. This difference is a result of how the 
spin couplings and energy denominators of the meson- 
exchange propagators are implemented in the kernels of 
the two frameworks. 

Before ending this section, a final comment concern- 
ing the accuracy of the numerical methods employed is 
provided. The accuracy is measured using the unitarity 
condition (or optical theorem), 



T -T* =T\G - G f )T, 



(3-11) 



which is derived from the fact t hat V in Eq. ( |2.1g| ) is Her- 
mitian. Evaluating Eq. Q3.ll ) between two-pion states, 
one obtains an equation that relates the two-pion flux 
missing from the forward direction to the one-, two-, and 
three-body outgoing flux observed leaving the scattering 
center. This relation provides a sensitive check of the 
numerical methods employed. 

The fraction of lost two-pion flux observed as outgoing 
two- or three-particle states is shown versus the expo- 
nent a of the adiabatic scale e = 10~ a GcV in Fig. |lj 
for CM energy E — 1400 GeV. For any energy E greater 
than the two-pion scattering threshold, there can be no 
stable one-particle state. It follows that no one-particle 
outgoing flux can be observed. In Fig. HJ, the regions 
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FIG. 14. Fraction of incoming two-pion flux at energy 
E — 1400 MeV in the isoscalar-scalar channel, appearing as 
"two-body" or "three-body" final states, or lost to unitarity 
violations. The fractions are shown versus the exponent a of 
the adiabatic scale e = 10 _a GeV. 



labeled "two-body" and "three-body" represent the frac- 
tions of the missing two-pion flux that appear as outgoing 
two-body (7T7r or KK) and three-body (irirp) states, re- 
spectively. Since the Hilbert space is restricted to one-, 
two- and three-body states, this should account for all 
the flux. However, in practice, the numerical methods 
employed introduce violations to the unitarity condition 
of Eq. (3.11). The region below the solid curve in Fig. IT3 



represents the fraction of flux that completely disappears 
from the theory; such a loss of flux violates unitarity. One 
observes that for values of e < 10~ 7 GeV, the violation 
is less than 10~ 6 of the outgoing flux and is therefore 
negligible. In the application to tttt scattering described 
above, the value of e = 10^ 12 GeV (or a = 12) was em- 
ployed. 



IV. SUMMARY AND FUTURE DIRECTIONS 

Herein, a framework suitable for the description of non- 
perturbative hadron scattering, based on the Bakamjian- 
Thomas formulation of relativistic quantum mechanics 
is introduced. In Sec. ||, it is shown that by including 
the interactions into the free mass operator M, one can 
ensure that observables calculated in the framework are 
Lorentz coviariant. When the Hilbert space is truncated 
to contain only one-, two- and three-body states, the 
resulting Lippmann-Schwinger equations (LSEs) form a 
closed set of coupled integral equations. The solution of 
these integral equations is obtained numerically. 

A significant improvement of this framework over ear- 
lier work is that the full effect of three-body states have 
been included. The three-body Green functions appear 
only in the two-body scattering kernel K and the two- 
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body self-energy E, and for energies E above the three- 
body thresholds, unitary branch cuts associated with 
these thresholds appear both in K and S. It follows 
that the kernel K and two-body self-energy E are com- 
plex functions of the energy E. The appearance of such 
three-body branch cuts provides the means for important 
dynamical effects, such as three-hadron production and 
decays into three-hadron states, which are automatically 
accounted for in our framework. Such effects have hith- 
erto been ignored for the most part in previous studies. 

The inclusion of three-body cuts into the integral equa- 
tions requires appropriate numerical methods to be em- 
ployed. These methods, necessary to solve the coupled 
set of LSEs are shown to maintain the unitarity of the 
theory to better than one part in one million. To demon- 
strate the utility of the framework, a preliminary study of 
7T7T scattering is car ried o ut in Sec. III. The simple model, 
introduced in Sec. 
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is able to provide an excellent 
description of the rnr phase shifts and inelasticities. 

The main purpose for developing this framework is 
to provide a means of incorporating the dynamics of 
low-momentum transfer, final-state interactions into the 
study of hadronic processes for energies up to a couple of 
GeV. In this energy region, the comparison between ex- 
perimental data and theorectical predictions from models 
of QCD for quark and gluon dynamics, are often made 
difficult due to the presence of final-state interactions. 
The soft rescattering of final-state hadrons tends to mask 
the QCD dynamics of interest. The framework developed 
herein provides a tractable means to incorporate the ef- 
fects of final-state interactions into studies of hadronic 
phenomena. 

Towards this end, the framework is constructed to be 
a consistent extension for the constituent quark model. 
It provides a means to unquench the quark model by pro- 
viding for hadron loops, multi-particle thresholds and the 
unitarity branch cuts associated with these. The result is 
the generation of complex-valued scattering amplitudes. 
It is an extension of the quark model, in that the quark 
model may be used to provide the elementary couplings 
and form factors for the hadronic interactions in V . The 
framework uses this real potential V to generate the full 
scattering solution. 

Future applications of the framework will focus on 
the dynamics of nucleon resonances and exotic, hybrid 
mesons that are the subject of current and proposed ex- 
periments at TJNAF. In the baryon case, for energies up 
to about 2 GeV, there are some well-known and strik- 
ing examples in tvN scattering for which three-body ef- 
fects are crucial in understanding the experimental ob- 
servables. For example, in the L2/.2.7 = P11 channel, the 
tvN inelasticity arising from the three-body tvtvN state is 
very large Q. It is likely that a complete understand- 
ing of the P\\ tvN scattering channel and the mysterious 
iV*(1440) resonance requires the full implementation of 
three-body unitary cuts that this framework provides. 

The effect of three-body cuts in exotic partial waves 
can also be very important. Most theoretical studies of 



hybrid meson decays have so far ignored effects of final 
state interactions. The best candidate for the exotic me- 
son, the tv 1 (1600) Q was found in the piv decay channel, 
which is predicted to be suppressed with respect to other 
two body decay channels, in particular the b\Tv p7| , p8| |. 
This shift of strength from b\Tv to piv could be explained 
by mixing with the three-body, lotvtv intermediate state 
which is believed to have a strong coupling to these two- 
meson channels. 
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